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ABSTRACT 

We investigate the effects of non-Gaussianity in the primordial density field on the reion- 
ization history. We rely on a semi-analytic method to describe the processes acting on the 
intergalactic medium (IGM), relating the distribution of the ionizing sources to that of dark 
matter haloes. Extending previous work in the literature, we consider models in which the 
primordial non-Gaussianity is measured by the dimensionless non-linearity parameter /nl, 
using the constraints recently obtained from cosmic microwave background data. We predict 
the ionized fraction and the optical depth at different cosmological epochs assuming two dif- 
ferent kinds of non-Gaussianity, characterized by a scale-independent and a scale-dependent 
/nl and comparing the results to those for the standard Gaussian scenario. We find that a 
positive /nl enhances the formation of high-mass haloes at early epochs, when reionization 
begins, and, as a consequence, the IGM ionized fraction can grow by a factor up to 5 with 
respect to the corresponding Gaussian model. The increase of the filling factor has a small 
impact on the reionization optical depth and is of order ~ 10 per cent if a scale-dependent 
non-Gaussianity is assumed. Our predictions for non-Gaussian models are in agreement with 
the latest WMAP results within the error bars, but a higher precision is required to constrain 
the scale dependence of non-Gaussianity. 
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1 INTRODUCTION 

Reionization marks a crucial event in the history of the universe, 
when the first sources of ultra-violet (UV) radiation ionize the neu- 
tral Intergalactic Medium (IGM) and affect the subsequent forma- 
tion of the cosmic structures. When reionization ends, the small 
amount of left neutral hydrogen is responsible for the absorption 
lines that we observe today in the spectra of far objects. How- 
ever, the way in which this complex phenomenon occurs is still 
not well understood and the most recent observations paint it as a 
spatially inhomogeneous and not istantaneous process. While the 
Gunn & Peterson trough of the high-z QSO spectra suggests a late 
epoch of reionization at z ~ 6 {Fa n et al. 2001; Becker et alj200ll : 
IWhite et alj|200j iFan et ai]|2006h . the very recent analysis of the 
5-year WMAP data on the cosmic microwave background (CMB) 
polarization shows an IGM optical depth r ~ 0.084 which is in 
better agreement wit h an earlier reionization redshift, z ~ 10.8 
dKomatsu et alj2 008). On the other hand, a late reionization end at 
z ~ 6 is also probed by the IGM temperature measured at z < 4 
dHui & Haiman1l2003l) and by the lack of evolution in the luminos- 



ity function of Lyman-q g alaxies betwe en z = 5.7 and z = 6.5 
dMalhotra & Rhoads 2004 see, however. loTa et al.1 d2008l) for evi- 
dences of a decline at high z). Overall, the present situation regard- 
ing reio nization at redshift z ~ 6 as probed by QSO spectra is still 
unclear dBecker et al]|2007h . 

Many analytic , semi-analyt i c and n umerical model s (see , 



e.g. iGnedinl 120001: Iciardi et alj l2003b|: IWvithe & Loebl J200: 



Barkana & Loeb 2004; Haiman & Holder 2003; Madau et al. 200' 
I Wvithe & Cer]|2007l ; Ichoudhurv & Ferrar j|2007l ; llliev et alj|2007l: 
iRicotti et alj|2008f) have been proposed to describe this poorly un- 
derstood reionization process. They basically relate the statistical 
properties and morphology of the ionized regions to the hierarchi- 
cal growth of the ionizing sources, making more or less detailed as- 
sumptions to describe the ionization and recombination processes 
acting on the IGM. Since the first sources of UV background radi- 
ation appear in the firstly formed dark matter haloes, which corre- 
spond to the highest peaks of the primordial density field, the reion- 
ization process is expected to strongly depend on the main param- 
eters describing the cosmological model and the power spectrum 
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of primordial density fluctuations. For instance, the possible pres- 
ence of an evolving component of dark energy can imprint signa- 
tures in the resulting morphology of the ionized r egions and change 
the time-scales of t he whole process (see, e.g. iMaio et al.ll200q ; 
ICrociani et alj20()3) . Also the nature and the statistical distribution 
of the primordial matter fluctuations can influence the reionization 
history. Although the standard scenario for the origin of the struc- 
tures assumes that the primordial perturbations are adiabatic and 
have a (almost) Gaussian distribution, small deviations from pri- 
mordial Gaussianity affect the dark halo counts, in the rare-event 
tail, thus also in the high peaks of the density field which originated 
collapsed objects at high z. 

Aim of this work is to investigate the effects of some level 
of non-Gaussianity in the primordial density field on the reioniza- 
tion history. We will make use of analytical techniques to describe 
the processes in action on the IGM. In particular, we will extend 
previous works in which the considered non-Gaussian models have 
density fluctuations described by a ren ormalized \ 2 probability dis- 
tribution with v degrees of freedom ( Avelino & Liddle 200(3) or 
by a modified Poiss on distribution with a given expectation value 
A dChen et al]|2003l) . Here we will adopt a more convenient way 
to introduce primordial non-Gaussianity, which has now become 
standard in the literature, based on the parameter /nl (see the next 
section for its definition). In particular, we will assume two differ- 
ent kinds of non-Gaussianity, characterized by a scale-independent 
and a scale-dependent /nl parameter. 

The paper is organised as follows. In Section [2] we introduce 
the main characteristics of the cosmological models with primor- 
dial non-Gaussianity considered here. Section [3] reviews the main 
assumptions of the analytical model adopted to describe the cosmic 
reionization process. The main results of our analysis are presented 
and discussed in Section[4] Finally in Section[5]we draw our con- 
clusions. 



2 MODELING PRIMORDIAL NON-GAUSSIANITY 

The main purpose of this work is to study the process of reioniza- 
tion under the assumption that the formation of the first ionizing 
sources is driven by the spherical collapse of overdense regions in 
a non-Gaussian primordial density field. The predicted reionization 
history will be compared to that obtained assuming the 'standard' 
model, with a Gaussian distribution of primordial perturbations, 
that will represent our 'reference' case. 

All the models considered in this work share the cosmological 
para meters suggeste d by the recent analysis of the 5-year WMAP 
data dKomatsuet al. 2008): a ACDM cosmology where the contri- 
butions to the present density parameter from dark matter, cosmo- 
logical constant and baryons are fi TO o = 0.279, f^Ao = 0.7214, 
f2;,o = 0.0461, respectively; the Hubble constant (in units of 100 
km/s/Mpc) is h — 0.701. The normalization of the cold dark matter 
power spectrum is fixed by assuming as = 0.817 and the primor- 
dial spectral index is taken to be n = 0.96. 

We will describe the level of primordial non-Gaussianity us- 
ing the dimensionless parameter /nl which weighs the quadratic 
correction to the linear Gaussian term in Bardeen's gauge-invariant 
potential $: 



$(x) = $ G (x) + /nl * (4(x) - <*g(x)» , 



(1) 



where $g (x) is a Gaussian random field and * denotes a convolu- 
tion. On scales smaller than the Hubble radius, $ is minus the usual 
Newtonian gravitational potential. With our convention, a positive 



value for /nl leads to a positive skewness for the distribution of 
the matter density fluctuations. 

As shown by eq.Q}, in general the non-Gaussian contribution 
to the gravitational potential $ can be written as a convolution be- 
tween a space- and/or shape-dependent /nl (x) and the quadratic 
term 4> G (x). The possible dependences of /nl are often neglected 
in the literature: this is done mostly for sake of simplificity, but it 
can motivated by the small r.m.s. value of In this case the bis- 
pectrum of the gravitational potential, defined as 

($(k 1 )$(k 2 )$(k 3 )> = (27r) 3 5 3 (k 1 ,k 2 ,k 3 )F s (fc 1 ,fc2,fc3) , 

(2) 

8 Z being Dirac's delta function, assumes a dependence on the 
wavenumbers called local shape, for which the term (fci , fc 2 , fc 3 ) 
can be expressed as 

F s (ki,k 2 ,k 3 )=2f^ h [P(ki)P(k2)+P(ki)P(k 3 )+P{k2)P(k 3 )] 

(3) 

In the previous equation P(k) = A$fc~ 3+ ' n_1 - ) represents the 
normalized power- spectrum of $. Bispectra which can be ex- 
pressed like in eq.10 are typical for models where the non- 
Gaussianity is produced outside the horizon or when the infiaton 
has a varying decay rate. It can be shown that eq.0 assumes the 
largest values for squeezed configurations, i.e. when one wavenum- 
ber is much smaller than the other two. 

Alternative models for primordial non-Gaussianity, based on 
a single field with higher derivative terms, predict a different shape 
for the bispectrum, having the so-called equilater shape. Its expres- 
sion can be still obtained using eq .Q, but replacing -F s (fci , fc 2 , fc 3 ) 
as follows dCreminelli et alj|20071) : 



F s (k 1 ,k 2 , k 3 ) = 6/ N L(fci, fc 2 , fc3)A $ 
2(fc!fc 2 fc3) 2(n - 1)/3 



{kik 2 



(kik 2 



+ 2 perm 



(fcife 2 fe 3 ) 2 



+ 



(kl /3 k 2 2 /3 k 3 )^ 



fcl kr~) fCq 



+ 



+ 5 perm 



(4) 



where we explicitly write the possible scale-dependence of /nl- 
As its name suggests, in this case the maximum amplitude of the 
bispectrum is where the wavenumbers are all equal. We refer to 
iLo Verde et al] (2008) for an extended discussion about the role 
of bispectrum s hapes in the parame trization of primordial non- 
Gaussianity and Bar tolo et alj d2004h for a review on the predic- 
tions for /nl in different inflationary models. 

At present, the strongest constraints on the parameter /nl 
are based on CMB data. Analysin g the 5-year temperatu re maps 
obtained by the satellite WMAP, iKomatsu et al . (2008) derived 
—9 < /nl < HI when the local shape for the bispectrum is 
assumed, and —151 < /nl < 253 for the equilateral shape. 
Both limits have been estimated at the 95 per cent confidence 
level. The fact that the presence of some amount of primordial 
non-Gaussianity alters the growth of density fluctuations and then 
the formation and evolution of cosmic structures, suggests that the 
large-scale structure (LSS) of the universe can be an alternative 
powerful probe for /nl, which has also the important advantage 
of being based mostly on three-dimensional data. Many theoretical 
studies, based both on analytic and numerical analyses, have inves- 
tigated the constraining capability of different observables like the 
abundances of virialised objects like cluste rs (Messina et al ]ll990l: 
| Moscardini et alj 1 19911: [Weinberg & Cold fl992l: iMatarrese et alj 

ogTE 



2000; IVerde et alj l200ol ; iMathis et alj |2004 iKang et alj 120071: 



The effects of primordial non-Gaussianity on the cosmological reionization 3 



Grossi et al]|2007l:lLo Verde et alj|2008l), halo bia sing jDalal et al.1 
20081 : Matarrese & Verdeil2008l;lMcDonaldl2008h . galaxy bispec- 
trum (Sefus atti & Komatsull2007l) , density mass field distribution 
dGrossi et alj|2008l) and topology dMatsubarall2003l: iHikage et al.l 
20081). integra t ed Sa chs-Wolfe effect dAfshordi & Tolled 1 120081 ; 
Carbo neet al.l 12008 ), low density intergalactic medium and 
the Lyman-q flux i Viel et al. 2008), 21 -centimeter fluctuations 
dCooray|2005l : |Pillepich et alj2007l) . and reionization, as discussed 
in this paper. In general, the application of these theoretical re- 
sults to real LSS data provided weaker constraints on /nl with 
respect to the CMB. The only e xception is the very recent anal- 
ysis made by ISlosar et al. (2008), who applied the bias formal- 
ism to a compendium of large-scale data, including the spectro- 
scopic and photometric luminous red galaxies from the Sloan Dig- 
ital Sky Survey (SDSS), the SDSS photometric quasars and the 
cross-correlation between galaxies and dark matter via Integrated 
Sachs-Wolfe effect. Considering the local shape only, they found 
—29 < /nl < 70 (at 95 per cent confidence level). It is impor- 
tant to notice that the scales probed by CMB and LSS are generally 
different and can give complementary information on /nl if the 
primordial non- Gaussianity is assume d to be scale-dependent (see 
the discussion in lLo Verde et al.| [2008). 

In this work, we consider non-Gaussian model with bispec- 
trum having both the local and equilateral shapes. We will use val- 
ues for /nl in the ran ge constrained by the 5-year WMAP results 
dKomatsu et ai]|2008l) , i.e. -9 < / NL < 111 and -151 < /nl < 
253 for local and equilateral shapes, respectively. In the last case, 
we also allow the non-Gauss ianity to vary w i th the scale, assuming 
the dependence proposed bv lLo Verde et alj d2008l) . namely: 



/nl(&1, &2, ks) — /NL 



fccMB 



(5) 



The normalisation of the previous relation is chosen in order to 
avoid violating the WMAP constraints: for this reason /nl rep- 
resents the equilateral parameter measured on the fccMB scale 
of O.O86/1/ Mpc, rough ly corresponding to largest multipole used 
by iKomatsu et al.l d2008h to estimate the non-Gaussianity in the 
WMAP data, I — 700. The slope a is a free parameter, assumed 
to be const ant, such that [q]_<g 1 between CMB and cluster scales. 
Following iLo Verde et all 12008) we consider small negative values 
for a, to enhance the non-Gaussianity on scales smaller than CMB. 
The resulting behaviour for the /nl parameter is shown in FigQ] 
where we assume a — 0, —0.1, —0.2, for the slope of the scale 
dependence, and /nl = —151, and /nl = 253 as pivoting values 
at the CMB scale, in agreement with the WMAP equilateral con- 
straints. It is evident from the plot that, with our assumption for the 
scale-dependence relation of non-Gaussianity, the absolute value 
of /nl at the scales relevant for the halo formation, and then for 
reionization, can be a factor 2-3 larger than the maximum amount 
directly derived from the CMB analysis: this can amplify the pos- 
sible effects of primordial non-Gaussianity. Furthermore, we note 
that the possible non-Gaussianity probes also are based on obser- 
vational data coming from different ranges of redshift. 

Deviations from Gaussianity influence the evolution of the 
density fluctuations, and affect the distribution of the virialized dark 
matter haloes at a given cosmological epoch. This reflects on the 
mass function: its high-mass tail is enhanced (reduced) in case of 
positive (negative) values of /nl- While in the Gaussian case, the 
IPress & Schechten d 19741) ( PS 74 hereafter) approach, together with 
its modern improvement s dLacev & Cold fl 993: Sheth & Tormenl 
Il999l : Ijenkins et al]|200ll : IWarren et alj|2006l) . represent powerful 



tools to model the evolution of the ionizing sources, for mildly 
non-Gaussian fields it is possible to m ake use of the analytic re- 
lation found by Matarrese et al. (2000) extending the PS74 formal- 
ism, which has been positively t ested against the re sults of high- 
resolution N-body simulations bv lGrossi et al.l d2007l) . If tig (M, z) 
represents the density of dark matter haloes with mass M at redshift 
z as obtained assuming Gaussi an initial conditions, he re modeled 
assuming the relation found bv lSheth & Tormenl (1999), the corre- 
sponding expression for a non-Gaussian models having the same 
cosmological parameters can be written as 



n NG (M,z) = F NG (M,z)n G (M,z) 
where the correction factor Fng is given by 



F NG (M, z) 



1 5l dS 3 (M) 
6 <5* d\na M 



+ 



exp 



Si S 3 



(6) 



(7) 



In the previous relation <5* = 8 c -\/l — S 3 (M)S c /3, 5 C represents 
the collapse threshold at z, a\i is the mass variance at z = and 
5*3 (M) is the normalized skewness of the primordial density field 
on mass scale M, namely Ss(M) — — f]s!L)J-'s{M)/a M . Then in 
order to compute the mass function for non-Gaussian models it is 
necessary to evaluate the third-order moment ^13, that depends on 
the bispectrum of the gravitational potential (<I > (ki) ( l>(k2) ( l > (k3)): 

„„ f dki f dk 2 f dk 3 



f dki r dk 2 r 

J (2tt)3 J (2^)3 J ' 



(27T) 3 

W(ki)W(k 2 )W(k 3 )F(k 1 )F(k2)F(h i ) x 
(S(ki)*(ka)*(ks)> 



(8) 



where W(k) is the Fourier transform of a spherical top-hat function 
on the mass scale M, F(k) = T(k)g(k), being T(k) the cold 
dark matter transfer function and g(k) = — 2(k/ Hq) 2 /(3Q m o) is 
required to go from the gravitational potential to the density via the 
Poisson equation. 

In Fig[2] we show, as a function of the halo mass, fi 3 and 
5*3, for models with both local and equilateral shapes. Both skew- 
ness parameters are given per unit non-Gaussianity parameter 
/nl; for the equilateral case we also consider the possibility of 
sca le-dependence f or the non-Gaussian term. As already shown 
by ILo Verde et alj d2008l) . the two classes of models give quite 
different predictions for both the amplitude and the mass depen- 
dence of the two considered quantities, but this discrepancy de- 
creases as the mass scale increases, since the local and the equilat- 
eral cases become more and more similar. The scale-dependence 
of /nl strongly affects the non-Gaussianity contribution for the 
smaller masses scales and this effect grows when higher negative a 
parameters are assumed. 

Inserting the values for S3 in eq.Q, we can estimate the ef- 
fects of primordial non-Gaussianity on the dark matter halo dis- 
tribution at the cosmological epochs relevant for the process of 
reionization. The results, shown in terms of ratio with respect to the 
Gaussian predictions, are shown in Fig(3] for both local and equilat- 
eral shapes (upper and lower panels, respectively). Here we adopt 
for the /nl parameter at the CMB scale the values corresponding to 
the 95 per cent confidence level, as derived from the 5-year WMAP 
data:/NL = —9, 111 for the local shape, and /nl = —151, 253 
for equilateral one. Since the mass density probability function is 
positively skewed in case of positive (/nl > 0) non-Gaussian 
contributions, the probability of overcoming the collapse thresh- 
old becomes higher. As a consequence, the formation of high-mass 
haloes is enhanced and anticipated when /nl > 0. Fig. [3] shows 
that the mass function can be increased by a factor of 10 at z = 13 
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Figure 1. The scale dependence of the non-Gaussianity parameter /nl for the models with equilateral shape considered in this paper. Results for two different 
choices of the pivoting value, /nl = —151 and /nl = 253, are shown by dashed and solid curves, respectively. Different values for the slope a have also 
been used: a = (black lines), a = —0.1 (red lines) and a = —0.2 (green lines). The shaded regions on the right show the scales probed by the SDSS (cyan) 
and the galaxy clusters (orange), while the blue region refers to the range excluded at 95 per cent confidence level by the CMB data (Komatsu et al. 20081). 
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Figure 2. The skewness [13 (top panel) and the normalized skewness S3 (bottom panel), given per unit /nl- The black dashed curve refers to the model 
with local shape, while the coloured solid curves present the results for the equilateral configuration with a = 0, —0.1, —0.2 (black, red and green lines, 
respectively). 
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for haloes with mass M ~ 10 Mq when compared to the stan- 
dard scenario. We should however remark that high-mass haloes 
(M > lO 9 Af0) at early cosmological epochs are rare events, as 
shown also by the small number density at z = 13 in the refer- 
ence case, n(> 1O 9 M ) < 5 x 1CT 3 /Mpc 3 . Then this effect is 
expected to have a little impact on integrated quantities as the total 
ionized fraction ot the IGM optical depth. Unlike the local model, 
the scale dependence of non-Gaussianity increases the abundance 
of the low-mass haloes by a factor of ~ 10 at z = 13 when com- 
pared to the standard case. The opposite ap plies for /nl < 0. 
As alr e ady noticed b y Matarrese et alj (2000) [see also lVerde et al.l 
d200lh : I Grossi et al. I 12QQ2I)], this effect is more evident at early 
cosmological epochs, exactly when the process of IGM ionization 
starts. For this reason a non-Gaussian distribution of the primor- 
dial density field can affect the way in which reionization occurs, 
leaving its imprints on it, as we will investigate in the next sections. 



3 AN ANALYTIC APPROACH TO COSMIC 
REIONIZATION 

In this section, we briefly review the main assumptions underly- 
ing the analytic model adopted to describe the process of cos- 
mi c reionization. Thi s mode l is based on the approach p r opose d 
by lAvelino & Liddlel d2006h [see also lHaiman & Holder! d2003l) : 
IChen et al] d2003l) for further details]; our implementation, how- 
ever, differs in some aspects which will be discussed later. 

In this model, the statistical properties of the ionized re- 
gions are related to the hierarchical growth of the ionizing sources 
through simple assumptions on how the galaxies ionize the IGM 
and on how the IGM recombines. A one-to-one correspondence 
between the distribution of galaxies and HII regions is established, 
such that a single galaxy of mass M ga i can ionize a region of mass 
Mhi = £M g ai. Here £ represents the ionization efficiency of the 
galaxy, and it is strictly dependent on the nature of the ionizing 
sources. We will take it as a constant, fixed in such a way that reion- 
ization ends at z = 6.5. 

Since at high z the cooling of the gas becomes efficient 
in haloes having a viri al temperature T 10 4 K, unlike in 
lAvelino & Liddlel J2006h . in our analysis we consider only type la 
(10 4 K ^ T ^ 9 x 10 4 K) and type lb (T > 9 x 10 4 K) haloes, 
neglecting the contribution of the type II sources, which would cor- 
respond to haloes with 400 K < T < 10 4 K. We recall that the 
distinction between the halo types is related to the way in which 
they impact the IGM: type la sources can grow only in neutral re- 
gions, while type lb haloes can appear also in ionized regions. Con- 
sequently they affect differently the ionization phases of IGM. 

The total collapsed fraction F 00 n [z) at different redshifts can 
be computed using eq.(6): 



F co \\{z) 



1 

Po 
1_ 

Po 



dM n NG (M, z) 



aO) 



dMn G {M,z)F NG (M,z) , 



(9) 



where po is the present-day matter density and M m i n (z) is the min- 
imum mass corresponding to the virial temperat ure T, which can 
be com puted by inverting the relation proposed by Bar kana & Loebl 
d200lh . namely: 



T = 1.98xl0 4 



10 



M 



10 8 Mb/i- 1 



2/3 



a. 



is?! 2 



1/3 

K. 

(10) 



In the previous equation, A c represents the virial overdensity at 
redshift z and is the matter density parameter at redshift z. 

Consequently, the collapsed fractions in la and lb haloes are 
given by 



Fcoll, Ib(z) — 
Fcoll, Ia(z) = 



1 f 00 

- / dMn NG (M,2) 

50 ^M min , /b (z) 



Po 
1 



dMn NG (M,z) - Feu, j& 0) 



(11) 



where M m i n jb and M m i n ,/ a are the minimum masses for lb and 
la sources, obtained using in eq. (10) T = 9 x 10 4 and 10 4 K, 
respectively. 

The action of the ionizing sources is smoothed down by the 
recombination of the IGM, here considered as a homogeneous gas. 
The recombination rate is linearly dependent on the IGM clump- 
ing factor Che =< w-h it > / < "-hi > 2 , for which, following 
lHaiman & Brvanl j200rj) . we assume a redshift evolution modeled 
as: 



C m {z) 



■9 



l + z 
7 



(12) 



being f3 a free parameter. As shown by lAvelino & Liddlel d2006l) . 
the predicted reionization history of the universe has significant un- 
certainties introduced by the poor knowledge of the z-dependence 
of the clumping factor, which cannot be robustly constrained even 
considering the 3-year WMAP results for the reionization optical 
depth. Since they found good consistency between predicted and 
observed optical depths irrespectively of the amount of primor- 
dial non-Gaussianity in the models they consider, we decide to set 
/3 = 0. In this case, the z-dependence of Che is neglected and 
Che = 10. The impact of the assumption of a constant clumping 
factor will be discussed later. 

The probability that a photon emitted at a given cosmological 
epoch Zi(ti) is still ionizing at z < z\ can be written as 



/>,/,/.. =oxp[ |- | 



(13) 



with ti = QBCHE«.Hi(ti)tf , being ob the recombination coeffi- 
cient of HI (= 2.6 x 10" 13 cm 3 /s at T = 10 4 K) and n m (z) = 
1.88f2 b0 /i 2 /0. 022(1 + z) 3 /cm 3 the hydrogen density at redshift z. 
Then, the filling factor Fhe(z) at a given cosmological epoch is 



Fhe(z) = / dz'C 



dF c( 



/ > J dFcoii, 



dz' 



-(z') + [1-Fhe(z')] x 



dz' 



{z')\P{z',z) 



(14) 



where the ionizing efficiency is assumed to be the same for the 
different types of haloes. We notice that the different nature of type 
lb and la sources appears in the right-hand side of eq.J14b, where 
the (1 — Fhi) factor explicitly considers that type la haloes form 
only in neutral regions. 

Finally, the HII filling factor allows us to estimate the reion- 
ization optical depth as follows: 



(z) = ott / dt' n c (t') 



1.08 cut 



dt_ 

d? 



dz' 1 



\ Y ) ea^L f hk( ,') . (i 5 ) 

4 / m D 
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Figure 3. The ratio between the dark matter halo mass functions for non-Gaussian and Gaussian models, computed at z = 6.5 (left panels), 2 = 9 
(central panels) and z = 13 (right panels). Top panels show the results for non-Gaussian models with local shape, where /m, = —9 (dashed lines) and 
/nl = HI (solid lines) are assumed. In the bottom panels, which refer to non-Gaussian models with equilateral shape, dashed and solid lines correspond to 
/ml = —151, 253, respectively; black, red and green lines refer to a = (i.e. no scale-dependence), a = —0.1 and a = —0.2 respectively. 



here <tt represents the cross-section of the Thompson scattering, 
n e is the free-electron density, c is the speed of light, Y is the He- 
lium mass fraction, ph(z) is the baryon density at redshift z and the 
factor 1.08 approximately accounts for the contribution of the Hel 
reionization, assuming that the HII and Hell fractions are equal and 
neglecting the effects of the Hell to Helll phase transition. 



4 RESULTS AND DISCUSSION 

In this section we present and discuss the main results of the appli- 
cation of the previous model under the assumption that the structure 
formation history is starting from a primordial non-Gaussian den- 
sity field. The corresponding predictions will be compared to the 
reionization scenario obtained for the Gaussian 'reference' case. 

First of all, using the minimum collapsed mass at each cos- 
mological epoch as de rived from the M-T relation proposed by 
Barka na & Loebl ( 2001), we can compute the total collapsed frac- 
tions for different kinds of sources, following eg. i ll It . The results 
are presented in Fig|4] where we show the collapsed fractions of 
la and lb haloes for the local and the equilateral cases here con- 
sidered and their redshift evolution compared to the predictions 
for the Gaussian scenario. We can notice that for both the local 
and the equilateral models, positive values for /nl produce an en- 
hancement of the collapsed fraction at a fixed redshift, given the 
larger probability for high-mass haloes in the primordial distribu- 
tion. When scale-dependent non-Gaussianity is considered, this ef- 
fect becomes higher as the power-law parameter a decreases. The 
opposite trend is observed when negative /nl values are assumed: 
in this case we obtain a reduction of the non-Gaussian collapsed 
fraction, due to the smaller probability for high-mass haloes at high 
z. 

The evolution of the collapsed fraction due to non-Gaussianity 



strongly affects the amount of ionized IGM at different cosmologi- 
cal epochs. As illustred in Fig|5] that shows the evolution of the fill- 
ing factor both for local and equilateral models, a positive primor- 
dial non-Gaussianity produces, if a scale-dependent /nl parameter 
is assumed, an increase of the ionized IGM density, that can be- 
come 5 times higher than that predicted in the Gaussian case. The 
IGM reionization results slower when negative values are consid- 
ered, since the smaller collapsed fraction produces a mild evolution 
of the filling factor, that is smaller compared to the Gaussian case, 
at every cosmological epoch. 

Finally in Fig[6]we show the effects of the primordial non- 
Gaussianity on the reionization optical depth r assuming both local 
and equilateral models. As a consequence of the evolution of the 
ionized fraction, in the initial phases of reionization the values of 
r for the models with positive /nl are higher than those predicted 
for the Gaussian case. In particular, when a scale-dependent non- 
Gaussianity is assumed, the change in r is higher than 10 per cent 
for a < —0.2. In this case, the non-Gaussian optical depth would 
be r ~ 0.083 at z ~ 30, that is still consistent with the last WMAP 
observations (r = 0.084 ± 0.016). Assuming a local shape with 
largest positive /nl here considered we find r ~ 0.079. Notice 
that for the Gaussian model we predict a reionization optical depth 
at z = 30 r = 0.078. 

We notice that our results are in qualitative agreement with 
the reionisation p icture resulting from the analysis performed by 
IChenet all (2003). However, a direct comparison cannot be done 
because of the differences in the considered cosmological parame- 
ters: in particular they assume the ones suggested by the first-year 
WMAP analysis, with a significantly higher power spectrum nor- 
malization (ag — 0.9). Furthermore, they parametrize the degree 
of non-Gaussianity using the expectation value of the considered 
modified Poisson distribution A, that is inversely proportional to 
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Figure 4. The collapsed fraction of la and lb ionizing sources. The top panels show the resulting F co \\ for the Gaussian reference case (dot-dashed line) and 
two extreme non-Gaussian models here considered, which adopt /nl = 111 with the local shape (solid line) and /nl = —151, 253 with a = —0.2 for the 
equilateral shape (green dotted and solid surves, respectively). The other panels show the ratio between the collapsed fractions of la (left panels) and lb (right 
panels) ionizing sources for non-Gaussian and Gaussian models as a function of redshift. In the middle panels the results for non-Gaussian models with local 
shape are shown, where /nl = —9 (dashed lines) and /nl = HI (solid lines) are assumed. In the bottom panels, which refer to non-Gaussian models with 
equilateral shape, dashed and solid lines correspond to /nl = —151, 253, respectively; black, red and green lines refer to a = (i.e. no scale-dependence), 
a = —0.1 and a = —0.2, respectively. 
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Figure 5. The left panel shows the evolution of the ionized fraction Fhi for the Gaussian and two different non-Gaussian models here considered, as in the top 
panels of Fig. [4] The two panels on the right show the ratio between the ionized fraction Fjji for non-Gaussian and Gaussian models as a function of redshift. 
In the left panel, the results for non-Gaussian models with local shape is shown, where /nl = —9 (dashed lines) and /nl = HI (solid lines) are assumed. 
In the right panel, which refers to non-Gaussian models with equilateral shape, dashed and solid lines correspond to /nl = —151, 253, respectively; black, 
red and green lines refer to a = (i.e. no scale-dependence), a = —0.1 and a = —0.2, respectively. 



the /nl parameter here used: A oc 1/ /nl • Their results show that a 
high amount of primordial non-Gaussianity (small A, i.e. high /nl) 
produces large deviations with respect to the Gaussian case in the 
ionized fraction and the IGM optical depth, in agreement with our 
results. We should, howev er, notice that for t he Gaussian reference 
model the estimate of r in lChen et al differs by 20 per cent 



with respect the value we predict using the 5-year WMAP cosmol- 
ogy- 

The predicted reionization optical depth (at z = 30) as a func- 
tion of the non-Gaussianity parameter /nl is shown in FigUJfor 
the local and equilateral cases. Indeed, the standard deviation on 
the current estimate of r from the last WMAP results (~ 20 per 
cent) does not allow us to strictly constrain the scale dependence of 
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Figure 6. The left panel shows the evolution of the IGM optical depth r for the Gaussian and two different non-Gaussian models here considered, as in the 
top panels of Fig. [4] The two panels on the right show the ratio between the reionization optical depth r for non-Gaussian and Gaussian models as a function 
of redshift. In the left panel the results for non-Gaussian models with local shape is shown, where /nl = — 9 (dashed lines) and /nl = 111 (solid lines) are 
assumed. In the right panel, which refers to non-Gaussian models with equilateral shape, dashed and solid lines correspond to /nl = —151, 253, respectively; 
black, red and green lines refer to a = (i.e. no scale-dependence), a = —0.1 and a = —0.2, respectively. 
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Figure 7. The reionization optical depth r at redshift z = 30 as predicted in different non-Gaussian scenarios is shown as a function of /nl ■ The dashed 
black curve refers to the local case, while the solid black, red and green lines correspond to the equilateral case with a = 0, a = —0.1 and a = —0.2, 
respectively. The blue and cyan rectangles represent the areas currently excluded by the WMAP observations for the local and equilateral cases, respectively. 



non-Gaussianity and, as shown auon &Riotto(2008), should 
affect the /nl value by a ~ 3 per cent and ~ 5 per cent uncer- 
tainty for the local and equilateral shapes, respectively. We estimate 
that in order to distinguish among the different scale dependences, 
a precision between ~ 1 per cent (/nl < 0) and ~ 8 per cent 
(/nl > 0) on the r measurement is required, using CMB-like ex- 
periments. In this sense, t he smaller standard deviati ons expected 
by Planck (~ 6 per cent) (Muk heriee & Liddlell2008h could better 
probe the reionization optical depth and would possibly constrain 
the non-Gaussianity parameter with an uncertainty of ~ 1 per cent 
and ~ 2 per cent in the local and equilaterals models, respectively. 
We should however keep in mind that the current WMAP estimate 
of r cited in this paper is affected by cosmic variance, since it is 
based on the bispectrum of the cosmic microwave background. Re- 
cent wor ks investigat e alternative methods to avoid this problem 
(see, e.g. lSeliakll2008l) . 

As a general warning, we stress that the description of the re- 
combination process in the analytic model here adopted relies on 
some simplifying assumptions which can affect the way in which 
reionization occurs. A deeper investigation of the IGM statistical 
distribution would require the analysis of suitable numerical sim- 



ulations having non-Gaussian initial conditions. As an example of 
possible biases introduced by the model uncertainties, in Fig[8]we 
show how a time-evolving clumping factor impacts the reioniza- 
tion optical depth. The res ults have been obta i ned by setting (3 = 2 
in eq . d 1 2b . as suggested bv lAvelino & Liddld (2006), for which the 
reionization optical depth is in agreement with the 3-year WMAP 
results for many degrees of non-Gaussianity. 

In order to estimate the relative importance of the details of 
the analytic modelling of the reionization process here adopted, we 
can compare the IGM optical depths predicted assuming (3 = 
(Fig[6j and (3 = 2 (Fig(8). As shown by the corresponding left 
panels, for the Gaussian model the effect on r due to the change 
in (3 and then to the unknown IGM physics is approximately of 
10 per cent. Viceversa, considering the non-Gaussian to Gaussian 
ratios (shown in the right panels of Figs|6]and[8j, changing (3 from 
to 2 produces negligible effects on the results. In fact we find 
that at high z the ratios obtained with (3 = and (3 = 2 differ 
by less than 1 per cent for the mildly non-Gaussian models and by 
~ 4 per cent for the non-Gaussian models with the most extreme 
scale-dependence (/nl = 253). By looking at the right panel of 
Fig. [8] we notice that for this last model the effect of introducing a 
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Figure 8. As in Fig[6] but assuming j3 = 2 in the redshift dependence of the clumping factor (see eq |12t . 



primordial non-Gaussianity can rise up to 10 per cent. This shows 
that the non-Gaussian to Gaussian ratios are only mildly affected 
by the choice of the recombination model: this justifies the fact that 
in this paper we preferred to show most of the results in terms of 
ratios in spite of absolute values. 

Finally, we remark that our approach, which is based on the 
PS74 formalism and its extensions, cannot fully account for the 
source clustering, that could indeed have relevant effects on the 
morphology of the HII regions. Since the source distribution de- 
pends on the clustering amplitude through the bias parameter, that 
is different if a primordia l non-Gaussianity i s considered (see, e.g . 
iMatarrese & Verdell2008l : iDalal et alj|2008l ; ICarbone et al] 2008). 
we expect that our results could be affected by this approximation, 
which can be improved only with suitable numerical simulations. 



5 CONCLUSIONS 

The aim of this work was to investigate how primordial non- 
Gaussianity may alter the reionization history when compared to 
the standard scenario based on Gaussian statistics. We have cho- 
sen a simple analytic method to describe the physical processes 
acting on the IGM to make predictions on the evolution of the 
ionized fraction and the reionization optical depth. Our work ex- 
tends previous analyses JChen et alj2003l ; lAvelino & Liddld2 006 ) 
based on simplified ways to introduce primordial non-Gaussianity, 
considering models motivated by inflation. In particular we as- 
sume two different hierarchical evolution scenarios for the ionizing 
sources, characterized by scale-independent and scale-dependent 
non-Gaussianity. All scenarios here considered are not violating the 
constraints coming from the recent analysis of the 5-year WMAP 
data. 

Our main conclusions can be summarized as follows: 

(1) non-Gaussianity affects the abundance of the dark matter 
haloes, since the formation of high-mass collapsed objects is en- 
hanced (reduced) when positive (negative) values for the /nl pa- 
rameter are assumed. This effect is more evident at earlier cosmo- 
logical epochs, exactly when reionization begins. 

(2) As a consequence, for positive primordial non-Gaussianity, 
the collapsed fraction in type la and lb sources is higher than for 
the Gaussian case at the same epoch, and the difference increases 
with z. The opposite result applies when /nl < is assumed. 

(3) The IGM filling factor is higher and its evolution is faster 
than in the Gaussian scenario if a positive /nl is assumed. This 
effect is enhanced for a scale-dependent non-Gaussianity, that can 



produce a 5 times higher Fjjm with respect to the Gaussian case, 
at early cosmological epochs. Viceversa the filling factor is smaller 
and has a mild redshift evolution for negative /nl • 

(4) Both local and equilateral non-Gaussianity have a small (less 
than 10 per cent) impact on the reionization optical depth, but the 
effect is enhanced assuming a scale-dependent /nl parameter. 

We finally remark that our predictions of the reionization op- 
tical depth in non-Gaussian cosmologies are in agreement with that 
estimated by 5-year WMAP analysis within la error bars and a 
precision higher than that of WMAP is required to constrain non- 
Gaussianity and its scale dependence. Ideally one would simulate 
reionization in non-Gaussian cosmological models using hydrody- 
namical simulations that incorporate all the relevant physical pro- 
cesses in a consistent framework (most importantly radiative trans- 
fer effects in the IGM). However, such approach is very time con- 
suming due to the large box size and the high resolution required 
to simulate large volumes and, at the same time, the physics of 
the sources of radiation. For this reason, some approximate semi- 
analytic schemes such as the one presentend here are still useful 
especially when calibrated on the more robust results of the hydro- 
dynamical runs. 
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